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1  Introduction 

Frequency  domain  substructure  synthesis  is  a  modeling  technique  that  enables  the 
prediction  of  a  combined  response  of  individual  structures  using  experimentally  measured  or 
numerically  predicted  frequency  response  functions  (FRFs).  The  traditional  synthesis  algorithm 
[1,2]  operates  on  component  impedances  and  thus  generally  requires  several  matrix  inversions. 
An  improved  algorithm,  developed  by  Jetmundsen  et  al.  [3],  requires  a  single  matrix  inversion 
with  a  completely  arbitrary  interface  definition  that  can  easily  incorporate  connection 
impedances.  The  main  limitations  of  the  method  are  the  large  data  requirements  and  sensitivity 
to  data  truncation. 

The  utility  of  this  technique  is  demonstrated  through  a  comparison  of  synthesized  and 
measured  admittances  of  an  edge-stiffened  plate  with  attached  equipment  (see  Figure  1.1).  The 
plate  mobilities  are  obtained  from  a  numerical  analysis  because  of  the  ability  to  accurately  model 
this  structure  using  a  finite  element  representation.  The  attachments  are  characterized 
experimentally  because  of  their  complexity.  The  following  sections  describe  the  synthesis 
technique  and  show  numerical  and  experimental  results  for  the  plate  and  equipment. 


Figure  1.1.  Edge-stiffened  plate  with  complex  attachments. 

2  Frequency  Domain  Impedance  Coupling  Technique 

The  frequency  domain  impedance  coupling  technique  represents  a  method  to  predict  the 
response  of  a  structure  using  substructure  dynamic  information.  Other  substructure  analysis 
techniques  are  available  and  are  commonly  used.  Of  these,  the  most  popular  are  finite  element 
substructuring  [1,4]  and  component  modes  synthesis  [1,2].  Finite  element  substructuring 
combines  structural  models  of  individual  structures  through  reduced  representations  of  the 
matrices,  while  component  modes  synthesis  combines  substructure  eigenvectors  to  compute 
eigenvectors  and  eigenvalues  of  the  combined  structure.  Often,  component  modes  synthesis  is 
used  in  a  finite  element  substructure  analysis  to  improve  the  accuracy  of  dynamic  predictions. 
Both  techniques  are  well  suited  for  numerical  analyses  where  the  mass  and  stiffness  matrices  can 
be  defined,  but  may  have  significant  shortcomings  when  dealing  with  experimental  data.  The 
frequency  domain  impedance  coupling  technique  can  be  used  for  numerical,  experimental,  or  a 
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combination  of  these  because  the  technique  is  based  on  component  mobilities  which  can  be 
directly  measured  or  numerically  computed. 

2.1  Traditional  FRF  Impedance  Coupling  Technique 

Impedance  coupling  via  frequency  response  functions  requires  knowledge  of  the  mobility 
matrices,  Y,  for  each  of  the  substructures  being  combined.  A  mobility  matrix  represents  the 
inverse  of  an  impedance  matrix  (Y  =  Z”1 ,  where  Z  is  the  structural  impedance),  which  is  defined 
as  follows: 

{p}=[z]-{x},  (1) 

where  {/?}  is  the  load  vector  and  {  x  }  is  the  velocity  vector. 

In  general,  the  mobility  matrices  are  6N  x  6N  complex-valued  matrices,  where  N  is  the 
number  of  nodes.  The  mobility  matrix  for  a  single  node  is: 
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where  the  coordinates  are  defined  in  Figure  2.1.  As  shown  in  Equation  2,  75%  of  the 
components  in  the  mobility  matrix  involve  a  rotational  coordinate,  an  important  detail  that  is 
discussed  later. 
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Figure  2.1.  Coordinate  system  used  for  the  definition  of  the  mobility  matrix. 


Mobility  matrices  similar  to  the  matrix  shown  in  Equation  2  can  be  defined  for  all  nodes 
of  each  structure  that  is  to  be  coupled.  Furthermore,  these  matrices  can  be  partitioned  according 
to  the  internal  (i)  and  coupling  (c)  DOFs.  The  partitioned  mobility  matrices  for  two 
substructures,  A  and  B,  for  frequency  co  are: 


2 


y»= 


(3) 


y‘  M 


- 1 

l - 

ya 

i  ya 

L  ci 

cc  J 

OQ 

1C 

i - 

yb 

cc  J 

(4) 


The  classical  derivation  of  the  FRF  substructure  synthesis  technique  employs  impedance 
matrices  instead  of  mobility  matrices  [1,2].  With  these  and  compatibility  of  displacements  and 
force  equilibrium  at  the  coupling  DOFs: 


u 


A 

C 


=  U 


B 

c 


fa+fb=f 

c  c  c 
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the  following  relation  for  the  coupling  of  substructures  A  and  B  can  be  derived  (see  Reference  2 
for  a  complete  derivation): 


Y  = 


0 


zA 

1C 

0 

ZA  +ZB 

ZB 

cc  cc 

Cl 

zB 

1C 

ZB 

ll 

(6) 


where  Y*  is  the  mobility  matrix  for  the  coupled  system.  As  seen  in  this  relation  (for  two 
structures),  such  a  coupling  process  requires  three  matrix  inversions:  one  to  obtain  the  impedance 
matrix  of  substructure  A  from  a  measured  or  predicted  mobility  matrix,  one  for  the  substructure 
B  mobility  matrix,  and  a  third  inversion  for  assembled  impedance  matrix  to  obtain  the  mobility. 
In  general,  the  number  of  matrix  inversions  is  one  greater  than  the  number  of  substructures  being 
combined.  The  size  of  these  matrices  are  ha  x  ha  for  substructure  A,  ns  x  ns  for  substructure  B, 
and  n  x  n  for  the  combined  structure,  where  n  =  ha  +  ns  -  nc •  It  is  also  apparent  from  Equation 
6  that  point  and  cross  mobilities  (or  impedances)  are  required  for  all  DOFs  and  all  frequencies 
involved  in  the  analysis.  It  is  this  fact  that  often  limits  the  usefulness  of  the  impedance  coupling 
technique. 


2.2  Generalized  FRF  Impedance  Coupling  Technique 

The  large  computation  requirements  of  the  traditional  impedance  coupling  technique 
(shown  in  Equation  6)  can  be  avoided  by  using  the  generalized  approach  developed  by 
Jetmundsen  et  al.  [3,5].  The  formulation  for  two  substructures,  A  and  B,  is  (see  Reference  5  for  a 
derivation): 
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It  is  apparent  from  Equation  7  that  only  a  single  matrix  inversion  of  size  nc  x  nc  is  required  for 
this  coupling  algorithm,  which  is  a  significant  improvement  over  the  traditional  approach.  The 
theory  has  been  generalized  to  incorporate  Boolean  mapping  matrices,  called  connectivity 
matrices,  that  define  the  structural  interconnections  [5].  The  generalized  relation  is: 
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(8) 


where  M  is  the  Boolean  mapping  matrix  defining  the  interfacial  sign  convention  for  the  coupling 
forces  and  their  reactions,  oris  a  substructure  identifier  in  addition  to  a  set  of  internal  DOFs,  yis 
an  interface  identifier  that  implies  an  interface  DOF  set,  ®  represents  a  matrix  element-by¬ 
element  multiplication,  and  N  is  the  number  of  substructures.  The  notation  is  best  understood 
with  the  use  of  an  example. 

Consider  the  structure  shown  in  Figure  2.2,  which  consists  of  three  substructures — A,  B, 
and  C — and  three  interfaces — I,  J,  and  K.  Each  column  of  the  mapping  matrix  corresponds  to  an 
interface,  while  each  row  corresponds  to  a  substructure.  A  positive  direction  is  arbitrarily 
assigned  to  each  interface,  as  shown  by  the  positive  and  negative  signs  adjacent  to  the  interface 
arrows  in  Figure  2.2.  The  arbitrary  sign  convention  is  applied  to  the  interfaces  for  the  purpose  of 
constructing  the  mapping  matrix.  For  each  substructure  (i.e.,  each  row  of  the  matrix),  a  “+1”  or 
“-1”  is  placed  in  the  columns  corresponding  to  an  interface  that  connects  to  the  substructure. 
The  decision  to  use  a  “+1”  or  “-1”  is  based  upon  the  chosen  sign  convention.  The  mapping 
matrix  corresponding  to  the  structure  in  Figure  2.2  is: 
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The  mobility  matrix,  Yaa,  is  defined  as: 
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Figure  2.2.  Frequency  domain  substructure  synthesis  diagram. 
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where  a  represents  the  internal  degrees  of  freedom  for  substructure  A,  b  are  the  internal  degrees 
of  freedom  for  substructure  B,  and  so  on.  These  DOF  sets  are  the  DOFs  for  which  mobility 
information  will  be  computed  in  the  synthesized  mobility  matrix.  If  response  information  is 
required  for  any  of  the  interface  degrees  of  freedom,  they  can  be  combined  in  these  DOF  sets. 
However,  each  interface  DOF  can  only  reside  in  one  set.  For  instance,  if  results  are  required  for 
the  interface  DOFs  associated  with  the  I  interface,  these  DOFs  can  be  included  in  either  the  a  or 
the  b  set,  but  not  both.  Note  that  the  order  of  the  DOF  sets  in  a  =  {a,  b,  c}  correspond  to  the 
order  of  the  substructures  in  the  rows  of  the  mapping  matrix  (Equation  9).  The  cross  mobility 
matrix  between  interface  and  internal  DOFs  is  defined  as: 
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Here,  the  i,  j,  and  k  DOF  sets  represent  the  DOFs  associated  with  the  7,  J,  and  K,  interfaces, 
respectively.  The  order  of  the  interface  DOF  sets,  y=  {i,j,  k},  corresponds  to  the  columns  of  the 
mapping  matrix.  When  multiplied  by  the  mapping  matrix,  M,  using  an  element  by  element 
multiplication,  some  of  the  components  in  this  matrix  are  nullified  (e.g.,  YAk  becomes  zero 
because  the  K  interface  does  not  involve  the  A  substructure): 
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The  final  part  of  Equation  8  that  requires  explanation  is  the  summation  term: 

<i3) 

9=1 

Here,  a  subscript  of  the  mapping  matrix  represents  a  row  of  the  matrix.  For  this  example,  M\  = 
[1  -1  0].  The  summation  from  q  =  1  to  TV  is  a  summation  over  all  substructures,  where  the  order 
corresponds  to  the  rows  of  the  mapping  matrix  (e.g.,  q  =  1  is  the  first  row  of  M  which 
corresponds  to  substructure  A,  q  =  2  is  the  second  row  of  M  which  corresponds  to  substructure  B, 
and  so  on).  With  this  definition,  we  have  Y'n  =  YA  where  ^is  the  internal  DOF  set  previously 
defined. 

While  the  theory  for  performing  substructure  synthesis  analyses  is  exact,  the  accuracy  is 
often  limited  by  incompleteness  or  inaccuracies  of  the  mobility  matrices.  These  deficiencies  can 
be  a  result  of  numerous  sources.  The  major  sources,  however,  are  often  a  result  of  mode 
truncation  or  lack  of  rotational  degrees  of  freedom  (RDOFs).  Both  of  these  sources  of  error  are 
discussed  next. 

2.3  Mode  Truncation  Errors 

Mode  truncation  errors  are  a  result  of  generating  frequency  response  functions  from  an 
incomplete  modal  database.  This  modal  database  can  be  derived  from  either  a  numerical  model 
or  experimental  measurements.  In  any  case,  the  errors  result  from  using  only  a  partial  set  of 
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eigenvectors  in  a  modal  summation  solution  for  a  system’s  response.  For  example,  a  finite 
element  model  may  consist  of  N  DOFs  from  which  N  mode  shapes  can  be  derived.  However, 
due  to  practical  limitations,  usually  only  the  first  M  modes  are  retained  (M  <  N).  Such  a 
truncation  generally  has  a  negligible  effect  on  the  response  at  resonance  frequencies,  but  does 
shift  the  antiresonances  upward  in  frequency.  (The  reason  for  the  upward  shift  in  antiresonance 
frequency  can  be  explained  by  considering  the  residual  term  that  is  discussed  later.)  While  this 
effect  may  have  a  negligible  impact  on  most  applications  of  the  FRFs,  substructure  synthesis 
methods  can  be  particularly  sensitive  to  it  [2,6]. 

The  effect  of  mode  truncation  is  shown  in  the  plots  of  Figure  2.3.  This  figure  shows 
results  from  modal  summation  solutions  for  a  cantilever  beam  finite  element  model  (structure  0 
in  Figure  2.4).  The  figure  indicates  the  increasing  errors  at  the  antiresonance  locations 
associated  with  the  drive  point  FRF  near  the  tip  of  the  beam  as  the  number  of  summed  modes 
decreases  from  7/  =  M  =  42  to  M  =  15  <  N.  The  consequence  of  these  errors  is  often  minimal  for 
applications  where  the  system  response  near  the  resonances  is  of  concern.  However,  for 
frequency  response  function  coupling,  errors  in  the  FRFs  at  all  frequencies  may  impact  the 
coupled  response.  This  is  demonstrated  by  synthesizing  the  response  of  the  beam  structure  of 
Figure  2.4  using  the  component  responses  computed  as  follows:  structure  0  response  computed 
using  a  direct  frequency  response  analysis  (equivalent  to  a  modal  summation  solution  using  all 
modes);  structure  0  response  using  a  modal  summation  with  modes  to  10  kHz;  structure  0 
response  computed  using  a  modal  summation  with  modes  10  kHz  and  a  residual  vector  to 
account  for  the  truncated  modes  (residual  vectors  are  discussed  below).  A  direct  frequency 
response  analysis  was  used  in  all  cases  for  structures  1  and  2.  The  drive  point  FRFs  at  grid  7  for 
the  combined  structure  are  shown  in  Figure  2.5.  This  figure  shows  excellent  agreement  between 
the  synthesized  responses  for  all  cases  except  the  case  using  a  modal  summation  and  no  residual 
vector,  wherein  differences  exist  in  both  the  resonance  and  antiresonance  frequencies.  Based  on 
these  results,  it  is  important  for  synthesized  results  that  errors  associated  with  modal  truncation 
be  minimized  whenever  possible. 


Cantilever  Beam  Model:  Modal  Summation  Solution 
(Drive  Point  Admittance  for  Grid  7,  ^-Direction) 


Figure  2.3.  Modal  summation  study  showing  mode  truncation  effects. 


6 


14 


Figure  2.4.  Beam  structure  example  problem. 


Cantilever  Beam  Model:  Frequency  Domain  Synthesis 
(Drive  Point  Admittance  for  Grid  7,  z- Direction) 


Figure  2.5.  Synthesis  example  showing  effect  of  mode  truncation. 

Because  it  is  often  not  practical  to  use  all  possible  mode  shapes  in  such  a  solution  due  to 
memory  limitations  for  numerical  modes  or  measurement  limitations  for  experimental  modes, 
the  effect  of  the  truncated  modes  can  be  approximated  with  the  use  of  residuals  [2].  Both  ‘low- 
frequency’  and  ‘high-frequency’  mode  truncation  can  be  accounted  for  using  residuals.  The 
truncated  low-frequency  modes  exhibit  a  mass-like  behavior  while  the  high-frequency  modes 
exhibit  a  stiffness-like  behavior  [2].  If  the  modal  summation  is  represented  in  terms  of  the  modal 

parameters  (resonance  frequencies  ( cor ),  damping  (rjr),  and  modal  constants  (rA/k)  for  a  general 
mobility  between  points  j  and  k),  the  modal  summation  has  the  following  form: 
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where  M  Rk  is  the  residual  mass,  K*k  is  the  residual  stiffness,  and  m\  and  m\  are  the  summation 
limits,  where  m\  indicates  truncation  of  low  order  modes  and  m2  -  mi  +  1  is  the  total  number  of 
modes  used  in  the  summation.  M*k  and  KR]k  are  dependent  on  the  applied  loading,  and  methods 

to  estimate  their  values  may  be  found  in  Reference  2.  Commercial  finite  element  solvers  also 
may  offer  methods  to  account  for  the  residual  terms  when  performing  a  modal  summation 
solution.  MSC/NASTRAN,  for  example,  can  compute  a  residual  vector  for  each  load  condition 
to  account  for  the  truncated  modes  by  including  the  bulk  data  entry  “PARAM,RESVEC,YES” 
[7].  The  process  involves  a  static  solution  of  the  model  to  obtain  deflection  shapes  that  are  then 
modified  or  discarded  to  make  them  linearly  independent  of  and  orthogonal  to  the  mode  shapes. 
The  residual  vectors  are  computed  and  appended  to  the  mode  shape  array  for  use  in  subsequent 
modal  summation  solutions. 

2.4  Errors  Caused  by  Rotational  Degree  of  Freedom  Truncation 

The  second  major  source  of  error  in  frequency  response  function  coupling  is  associated 
with  neglecting  components  of  the  mobility  matrices.  Often  the  components  not  included  are  the 
mobilities  associated  with  rotational  coordinates  [8].  As  discussed  previously,  75%  of  the 
complete  mobility  matrix  involves  a  rotational  coordinate.  The  importance  of  these  mobilities 
depends  on  the  components  being  coupled  and  the  interface  DOF  placement.  For  example,  a 
structure  with  a  broad  interface  that  is  connected  using  several  DOFs,  such  as  the  structure 
shown  in  Figure  2.6a,  may  not  require  the  use  of  the  rotational  terms  because  the  translational 
terms  inherently  provide  the  rotational  information  [1],  However,  a  structure  that  has  point  or 
collinear  interface  DOFs  (see  Figure  2.6b)  may  require  all  rotational  terms  or  rotational  terms 
about  the  collinear  axis,  respectively. 


Figure  2.6.  Choices  of  interface  DOFs  to  mitigate  effects  of  not  including  RDOF 

mobilities:  (a)  interface  DOFs  may  adequately  capture  rotational  effects,  and 
(b)  collinear  DOF  placement  may  require  RDOF  mobilities. 

The  need  for  rotational  degrees  of  freedom  is  demonstrated  using  the  plate/mass  model 
shown  in  Figure  2.7.  Here,  the  mass  (structure  B)  is  joined  to  the  plate  (structure  A)  at  two 
collinear  points  on  the  interface.  The  combined  response  is  computed  in  three  ways:  1)  Nastran 
solution  to  the  combined  model  (i.e.,  the  true  solution),  2)  frequency  domain  synthesis  using 
translational  DOFs  (x,  y,  and  z )  only,  and  3)  frequency  domain  synthesis  using  DOFs  x,  y,  z,  and 
(3.  The  drive  point  response  at  one  of  the  interface  grids  for  each  of  these  cases  is  shown  in 
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Figure  2.8.  This  figure  shows  good  agreement  between  cases  1  and  3  for  both  resonance 
frequencies  and  amplitudes  across  the  entire  spectrum.  However,  the  results  for  case  2  indicate  a 
missed  resonance  at  155.  Hz,  incorrect  resonance  and  antiresonance  frequencies  at  several 
locations  (e.g.,  410.  Hz),  and  different  amplitudes  at  several  frequencies. 


Figure  2.7.  Plate  example  for  rotational  DOF  truncation. 


Plate  and  Mass  Model:  Effect  of  Rotational  DOFs  on  Synthesis 
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Figure  2.8  Results  for  plate  model  showing  effect  of  RDOF. 

Rotational  DOF  mobilities  are  usually  neglected  due  to  the  difficulty  associated  with 
measuring  rotational  moments  and  accelerations  at  a  point.  While  various  attempts  have  been 
made  to  develop  an  accurate  method  for  measuring  these  mobilities,  the  most  commonly  used 
approach  is  to  infer  the  rotational  information  from  measured  translational  information  using  a 
finite  difference  approach  [9].  ANSI  S2. 34-1984  defines  standard  methods  for  measuring  or 
estimating  the  RDOF  mobilities  using  either  translational  measurements  or  rotational 
measurements  involving  specialized  fixtures  [10]. 
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Figure  2.9.  Translational  mobility  measurements  for  the  Finite  difference  RDOF 
approximation. 

The  finite  difference  approach  involves  the  measurement  of  translational  point  and  cross 
mobilities  as  shown  in  Figure  2.9.  The  number  of  required  measurements  is  governed  by  the 
desired  accuracy  of  the  results.  The  order  of  accuracy1  is  one  less  than  the  number  of 
measurement  points  when  the  traditional  finite  difference  formulae  are  employed.  For  example, 
first  order  accuracy  requires  the  measurement  of  mobilities  for  two  locations,  second  order 
requires  three  locations,  and  so  on.  A  compact  implementation  of  this  method  has  been 
implemented  by  Duarte  and  Ewins  [9].  Their  method  employs  the  following  relation: 
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where  Hineas  is  the  measured  mobility  matrix,  Hest  is  the  estimated  mobility  matrix  containing 
RDOF  terms,  and  T  is  the  transformation  matrix  that  is  defined  according  to  the  approximation 
order.  The  transformation  matrices  for  first  and  second  order  schemes  are  shown  in  Table  2.1. 
These  matrices  require  a  constant  spacing  of  the  measurement  locations  similar  to  the  spacing 
shown  in  Figure  2.9.  The  forward  and  backward  differencing  schemes  shown  in  this  table  are 
required  for  estimating  mobility  information  near  the  edges  of  a  structure.  Transformation 

Table  2.1.  First  and  Second  order  finite  difference  transformation  matrices. 


First  Order 

Second  Order 

Forward 

'  0  1 

1  Is  -l/s_ 

1 

2s 

"0  0  2s' 

-1  4  -3 

Central 

N/A 

1 

2s 

'  0  25  O' 

-l  0  1 

Backward 

1 - 1 

1 

o 

1 _ 1 

1 

2s 

'0  0  2s' 

1  -4  3  _ 

1  The  prediction  error  is  related  to  the  sample  spacing  as  follows:  Vh  order  has  accuracy  on  the  order  of  5"v  ,  where  S 
is  the  sample  spacing. 
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matrices  for  unequal  spacing  and  approximation  orders  other  than  first  or  second  can  be 
developed  using  the  general  approach  described  in  Reference  10.  Here,  the  RDOFs  and 
translational  mobilities  are  related  in  terms  of  the  partial  derivatives  of  the  translational  mobility 
at  the  point  of  interest.  Formulae  to  approximate  the  partial  derivatives  are  then  developed  based 
on  the  finite  difference  approach. 

3  Synthesis  of  Plate  and  Attachment  Response 

As  stated  previously,  the  frequency  domain  substructure  synthesis  technique  is  aptly 
suited  to  accommodate  experimental,  numerical,  or  any  combination  of  these  data  sources.  The 
edge-stiffened  plate  (see  Figure  3.1)  with  attached  electronic  equipment  (DC  fan  and  switching 
power  supply  shown  in  Figure  3.2)  represents  a  case  where  a  combination  of  numerical  and 
experimental  data  is  most  appropriate.  The  advantage  of  using  a  numerical  model  over  an 
experimental  characterization  is  mainly  due  to  the  large  data  requirements  of  the  synthesis 
process.  Numerically,  the  point  and  cross  mobilities  are  easily  acquired  through  multiple  load 
cases.  Acquiring  similar  data  experimentally  could  require  considerable  effort.  However, 
depending  on  the  component  complexity  it  may  be  more  effective  to  experimentally  characterize 
a  component  rather  than  invest  a  significant  amount  of  time  developing  an  accurate  model  of  the 


Figure  3.2.  Attached  equipment  showing  the  mobility  measurement  locations:  (a)  the  DC 
Fan  and  (b)  the  switching  power  supply. 
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structure.  The  DC  fan  and  switching  power  supply  are  two  such  components  that  would  require 
significantly  more  time  to  model  than  to  experimentally  characterize.  Based  on  this,  the  plate 
mobilities  are  determined  from  a  FE  model  while  the  attachments  are  characterized 
experimentally.  A  description  of  the  plate,  the  finite  element  model  of  the  plate,  and  the 
electronic  equipment  is  provided  in  Reference  11. 

Even  though  the  plate  has  a  relatively  simple  construction,  the  predicted  mobility  matrix 
is  slightly  different  from  the  measured  mobility  matrix  (mainly  due  to  inaccuracies  in  the 
plate/stiffener  interface  representation  in  the  model).  The  main  differences  involve  the 
resonance  frequencies  and  damping  levels,  as  shown  by  differences  in  frequency  and  magnitude 
of  the  peaks  in  the  mean  surface  admittance  (see  Figure  3.3).  Because  the  performance  of  the 
synthesis  technique  is  being  investigated  here,  it  is  important  to  tune  the  plate  model  to  match  the 
experimental  results  such  that  the  differences  associated  with  modeling  errors  are  minimized. 
The  model  was  tuned  by  modifying  the  computed  resonance  frequencies  and  damping  to  match 
those  determined  experimentally.  This  process  involves  two  steps:  1)  an  identification  step 
where  the  numerical  and  experimental  modes  are  matched  based  on  a  modal  assurance  criteria 
(MAC)  analysis  (see  Reference  11)  between  the  numerical  and  experimental  results,  and  2)  an 
adjustment  step  where  the  numerical  resonance  frequencies  and  modal  damping  levels  are 
adjusted  to  match  the  experimental  results.  The  result  of  this  process  is  a  model  that  has  a 
resonant  character  that  accurately  represents  the  dynamic  characteristics  of  the  plate  up  to  a 
frequency  of  500  Hz.  The  mean  admittance  for  the  experimental  and  numerical  results  due  to  a 
common  drive  point  in  the  normal  direction  is  shown  in  Figure  3.4.  As  indicated  by  this  figure, 


Measured  and  FE  Results  for  Bare  Plate 


Figure  3.3.  Mean  measured  and  predicted  admittances  for  the  bare,  edge-stiffened  plate. 
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the  modified  numerical  and  experimental  data  sets  agree  well  in  terms  of  resonance  frequency 
location  and  peak  amplitude.  Note  that  the  resonances  below  ~50  Hz  are  rigid  body  modes  that 
result  from  the  mass  of  the  plate  and  the  stiffness  of  the  bubble  wrap  used  during  testing  (see 
Figure  3.1). 

The  next  step  in  the  analysis  is  to  define  the  interface  degrees  of  freedom.  The  narrow 
attachment  surfaces  of  the  equipment  prohibited  an  accurate  estimation  of  rotational  degrees  of 
freedom  at  the  attachment  points.  Therefore,  interface  points  were  located  at  the  four  comers  of 
each  attachment  surface  to  inherently  capture  the  rotational  aspects  of  the  structure  as  previously 
discussed.  In  addition,  only  the  translational  DOFs  normal  to  the  plate  surface  were  used  in  the 
analysis  because  of  the  difficulty  associated  with  measuring  the  in-plane  mobilities.  Because  the 
majority  of  the  plate  modes  involve  out-of-plane  motion,  excluding  the  in-plane  DOFs  should 
not  have  a  significant  impact  on  the  results.  The  attachment  locations  are  shown  in  Figure  3.2. 

The  mobility  matrices  for  the  bare  plate  were  obtained  by  performing  a  modal  frequency 
response  using  the  modified  eigenvalues  and  damping  levels  as  discussed  above.  Effects  of 
mode  truncation  were  minimized  by  using  eigenvectors  with  frequencies  up  to  4  kHz  and 
MSC/Nastran’s  residual  vector  method.  Point  and  cross  mobilities  for  all  of  the  attachment 
DOFs  were  obtained  through  the  use  of  several  load  cases. 


Measured  and  (Modified)  FE  Results  for  Bare  Plate 


Figure  3.4.  Mean  measured  and  modified  numerical  admittances  for  the  bare,  edge 
stiffened  plate. 
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The  mobility  matrices  for  the  attached  equipment  were  obtained  experimentally  using  a 
tap  test  (an  impact  hammer  to  measure  the  applied  force  and  accelerometers  to  measure  the 
response).  Results  from  these  measurements  for  the  translational  DOFs  are  shown  in  Figure  3.5 
for  the  DC  fan  and  in  Figure  3.6  for  the  switching  power  supply.  The  plots  in  these  figures  show 
the  highly  damped  equipment  mobilities.  The  impedance  coupling  technique  is  a  valid  candidate 
for  this  type  of  data  since  an  accurate  modal  model  (which  is  required  for  other  synthesis 
techniques)  would  be  difficult  to  derive.  It  should  be  noted  that  the  low  frequency  resonances 
(approximately  less  than  50  Hz)  shown  in  these  figures  are  rigid  body  modes  that  result  from 
resting  the  structures  on  bubble  wrap  during  the  testing. 

The  mobilities  of  the  combined  structure  are  predicted  using  the  measured  equipment 
mobilities  and  the  numerically  generated  mobilities  of  the  bare  plate.  The  mobilities  of  the  plate 
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Figure  3.5.  DC  fan  mobility  matrix  components  (referenced  to  the  grid  numbers  shown  in 
Figure  3.2a). 
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and  attachments  are  also  obtained  experimentally.  The  attachments  are  joined  to  the  plate  using 
J-B®  WELD  Kwik  Weld  adhesive  and  small  flat  washers  at  the  four  mounting  locations  of  each 
component  to  provide  point-like  attachments.  The  results  of  both  data  sets  are  provided  in 
Figure  3.7.  These  results  show  a  similar  shift  in  both  resonance  frequencies  and  amplitudes 
between  the  bare  plate  and  plate  with  equipment  mobilities.  Note  that  the  synthesized  results 
ignore  the  effect  of  rotational  mobilities  and  translation  for  the  in-plane  directions. 


200  400 


200  400 


200  400 


200  400 


200  400 


200  400 


200  400 

Frequency  (Hz) 


200  400 


Figure  3.6.  Switching  power  supply  mobility  matrix  components  (referenced  to  the  grid 
numbers  shown  in  Figure  3.2b). 
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Measured,  Modified  FE,  and  Synthesized  Results 


Figure  3.7.  Measured  and  predicted  bare  plate  and  combined  plate  and  equipment 
admittances. 

4  Summary  and  Conclusions 

The  generalized  frequency  domain  substructure  synthesis  technique  proposed  by 
Jetmundsen  et  al.  has  been  used  to  compute  the  combined  response  of  a  bare,  edge-stiffened 
plate  and  two  electronic  components.  The  FRF  coupling  technique  was  chosen  over  other 
coupling  methods  because  it  is  better  suited  to  combine  data  sets  of  varying  origin.  Furthermore, 
Jetmundsen’ s  generalized  FRF  technique  has  been  used  because  it  offers  a  significant 
computational  improvement  over  the  traditional  frequency  domain  technique.  The  generalized 
method  is  also  amenable  to  interfacial  impedances  and  offers  completely  arbitrary  interface 
definitions  for  the  combination  of  multiple  substructures. 

The  major  shortcomings  of  the  frequency  domain  techniques  are  identified  as  modal 
truncation  and  incomplete  mobility  matrices.  The  modal  truncation  effects  are  only  important 
when  mobilities  are  acquired  from  modal  summation  solutions,  and  these  effects  can  be 
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minimized  through  the  inclusion  of  residual  terms.  Incomplete  mobility  matrices  are  a  result  of 
the  current  technological  limitations  for  measuring  rotational  mobilities  at  a  point.  However, 
methods  are  available  for  estimating  these  mobilities  when  they  are  deemed  important. 

The  synthesis  results  for  the  plate  and  equipment  show  similar  shifts  in  resonance 
frequency  locations  and  amplitudes  relative  to  the  experimental  results  for  the  combined 
structures.  Differences  in  the  results  have  been  attributed  to  limitations  of  the  experimental 
characterization  of  the  equipment.  The  large  data  requirements  of  this  method  may  be  a  limiting 
factor  when  trying  to  predict  the  coupled  response  of  entire  surfaces,  which  is  required  for 
radiated  sound  predictions.  For  components  with  small  surface  areas  or  applications  requiring 
only  a  few  FRFs  (such  as  loss  factor  predictions),  the  large  data  requirements  may  not  be  an 
issue. 
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